Coherent states superpositions in cavity quantum electrodynamics with trapped ions 
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We investigate how superpositions of motional coherent states naturally arise in the dynamics of a two-level 
trapped ion coupled to the quantized field inside a cavity. We extend our considerations including a more 
realistic set up where the cavity is not ideal and photons may leak through its mirrors. We found that a detection 
of a photon outside the cavity would leave the ion in a pure state. The statistics of the ionic state still keeps some 
interference effects that might be observed in the weak coupling regime. 



I. INTRODUCTION 

' There has been a great deal of interest in the coherent manipulation of simple quantum systems flSIUli mainly to the high 
, degree of control necessary for the implementation of quantum information processing tasks j^lUla]. In particular, the study 
ps| of trapped ions interacting with laser beams has attracted much attention due to the significant experimental advances in the 
I , generation of quantum states in such a system |7, 8]. The interaction of trapped ions with laser beams is well understood in 
^ ■ terms of a semiclassical model with the electromagnetic field being treated as a c-number, but new features mikght be revealed 
due to the field quantization. The entanglement between photons and ions is a remarkable consequence of that quantization and 
its potential applications heve been motivating the experimental work in cavity quantum electrodynamics with trapped particles 
10. For instance, there have been reported schemes for the generation of specific entangled states such as Greenberger-Horne- 
^ ! ZeiHnger (GHZ) states O as well as Bell states 

One of the reasons for interest in studying and experimentally coupling photons with material particles comes from the fact 
1^ that in order for the quantum information processing to be used in its full extent, one should be able to inter-convert stationary 

• and flying qubits and also faithfully transmit the flying qubits between given positions. Those two statements are part of what 
Cn| ', it are know as DiVincenzo's requirements for the physical implementation of quantum computation and information IT2I1 . The 

entanglement present in the system consisting of cavities and trapped ions may be useful in the propagation of information 
carried by photons between two distant locations 1131. 

It is not just entangled states involving either two level systems orFock states of the electromagnetic field that find applications 
in quantum information. Another interesting class of nonclassical states with high potential applications is the one formed by 
Qh, linear superpositions of coherent states. This class of states has been considered for quantum teleportation logic gates 

• implementation |16, 17] and tests of local realism fl^ . for instance. In this paper, we show that superpositions of motional 
^ coherent states may be generated in the framework of cavity electrodynamics with trapped ions by letting the system evolve 
2 in the resonant carrier dynamics and by performing a measurement of the internal state of the ion. We apply the formalism of 
qh' quantum jumps to study the non ideal case including damping in the cavity and show that the detection of photons outside the 
• • . cavity could be used to generate nonclassical states of the motion of the trapped ion. More precisely, we show that the statistics 

. ^ of the generated state keeps track of the coherence displayed on the oscillatory behavior of the phonon number distribution 
^> and the variation of its width from Poissonian to sub or super- Poissonian |19!|. Although that is not a deterministic protocol (it 
^ ' depends on the random event of the leaking of a photon from the cavity), it might be of interest because it could be implemented 
in current experimental systems. Experiments involving trapped ions and optical cavity fields have been performed only in the 
weak coupling regime in which the cavity damping is stronger than the ion-cavity coupling 1^. 



II. MODEL HAMILTONIAN 

In this work we consider a single two-level ion trapped in a Paul trap and placed inside an optical cavity. The cavity mode 
couples to the ionic internal degrees of freedom {|e),|(7}} and the system Hamiltonian is given by lEoll 

H = hva) d + hujb^b + h-^&z 

+hg{a+ + d-){b^ + b)cosri{a) +a), (1) 

where a) [a) denotes the creation (annihilation) operator of the center-of-mass vibrational motion of the ion (frequency v), w{b) 
is the creation (annihilation) operator of photons in the field mode (frequency lo), a operators are the usual Pauli matrices for the 
two internal levels of the ion, ujq is the atomic frequency transition, g is the ion-field coupling constant, and r\ ~ 27rao / A is the 
Lamb-Dicke parameter, being oq the amplitude of the harmonic motion and A the wavelength of the cavity field. 
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For our purposes here we may work in the Lamb-Dicke regime (77 ^ 1), i.e., the situation in which the spatial extend of 
the motion of the trapped ion is much smaller than the wavelengh of the cavity field. In this regime, we may perform an 
approximation that simplifies the original Hamiltonian Q as follows 



cos?7(d^ + a) ~ 1 



(2) 



If we tune the hght field so that it exactly matches the atomic transition, i.e., ujq — ui = (carrier transition), we obtain the 
interaction Hamiltonian in the Lamb-Dicke regime, which, after discarding rapidly oscillating terms reads 



Hi = hg 



1 - 



77^(1 + 2a^a) 



a+b). 



(3) 



The resulting Hamiltonian in equation is similar to the Jaynes-Cummings Hamiltonian but having an effective coupling 
constant which in our case depends on the excitation number of the ionic oscillator, m ~ a^a. Such a dependence on the 
intensity has already been demonstrated I21I1 to be related to the occurrence of super-revivals (revivals taking place at long 
times) of the atomic inversion. 



III. RESULTS 

A. Generation of superpositions of motional states 

We now consider that the system is initially prepared in a way that the ion is in its excited state |e) (internal level), the cavity 
in the vacuum state |0)c, and the vibrational motion in the coherent state \a)v, i.e. \ip{Q)) = \a)y |0)c|e). Under the Hamiltonian 
Q, the state |'0(O)) evolves to 

\^{t)) ^ cos{gt[l^r^^{l + 2a^a)/2])\a),\0)c\e) 

sin {gt [1 - 7/2(1 + 2ata)/2]) |a),|l),|g). 

(4) 

We still have to apply the functions of the operator a^a in the coherent state \a)u. This may be easily done by moving to the 
Fock basis and the result is given by 

\i{;{t)) = [cos(w^t) -ism{u}^t) |0)c|e) 

+ [cos(a;^t) -ism{u,,t) \<^+)v] |l)c|g), 

(5) 

where = g{l — /2) and we |$±)u are general superpositions of coherent states given by 

_ |ae'^)„±|ae-'^). 
\'^±)v = ^ , (6) 

where we defined the time dependent real phase (f) ~ V^gt- The state Q is an entangled state involving superpositions of 
motional coherent states of the trapped ion, its internal electronic states, and Fock states of the cavity field. It is noteworthy that 
for interaction times given by tk ~ kir, with k being an integer number, the state of the system reduces to 

m)) = |<f+)„|0),|e) + |$_},|l),|g}. (7) 

One could then obtain a disentangled motional state by performing a measurement on the internal state of the ion. The exper- 
imental discrimination between the two electronic levels may be done using the very efficient electron shelving method 
Depending on the measurement outcome, the collapsed motional state may be either or \^-)v One of the main interest- 

ing characteristics of those superposition states is that their statistics are strongly sensitive to the value of the phase cj). The trivial 
case takes place when = 0, what leads the distribution P,n = | {m\^±)y\'^ (phonon statistics) to be Poissonian. However, it is 
well known that there are domains in which it can be either sub or super-Poissonian. As pointed out in fl^ . when the statistics is 
super-Poissonian, the distribution P„i displays an oscillatory behavior, being this a direct consequence of interference in phase 
space. Such a behavior is analogous to the oscillatory photon statistics of highly squeezed states t23J . Although similar super- 
position states may also be generated using classical fields l24il . the possibility of entanglement with hght is a unique feature 
related to the quantum nature of the electromagnetic field. 
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FIG. 1: Ratio between the exact coupling constant and the approximate one in the Lamb-Dicke regime. The approximation is valid on the 
region R{rj, in) ~ 1 where r) and m are small enough. 



FIG. 2: Schematic experimental setup. The system consists of a single trapped ion placed in a lossy cavity having a decay rate k. The detector 
D continuously monitors this decay channel. 



The scheme proposed here relies on a not very demanding initial preparation of the system. It requires the initial field to be in 
the vacuum state |0)c, i.e., there is no need to prepare or inject a coherent field state into the cavity. Additionally, the vibrational 
motion of the ion has to be prepared in a coherent state \a)y, whose experimental realization for a ^Bc^ ion trapped in a RF 
(Paul) trap has been already reported 0]. Regarding the internal ionic states, they need to be prepared in the excited state which 
can be achieved by the application of laser pulses, for instance. 

We would like to point out that the linear dependence of the ion-field coupling constant on the operator a) a is crucial for 
the generation of the states |<i>±)u. Therefore, it is very important to be aware of the limits where the parameter 77 and the 
initial magnitude a may be varied and still having the approximation (|3 valid. This limit is set by keeping the product rf'a^a 
small enough, what allows us to neglect higher order terms in the cosine expansion. If the Lamb-Dicke approximation was not 
performed, it would be necessary to work with the full nonlinear coupling constant A = (m| cos ri{a) +a)\m) = /'^L^^{rf'). 
For convenient values of the product rf'a^a this coupling constant reduces to Xld = 1 — if {1 + 2m)/2, that is the coupling 
constant we have used so far (Lamb-Dicke regime). In figureQlwe show the ratio between the exact and the approximate coupling 
constants, R{ri, m) = X/Xld- We see that there are ranges of values of ?/ and m for which i? w 1. Under such circumstances, 
the Lamb-Dicke approximation is valid and the generation protocol proposed here is applicable. 



We are now interested in the more realistic setting where the cavity is not ideal and one could detect photons leaking through 
its mirrors. The set up we have in mind is depicted in Fig|2] We still have a two-level trapped ion interacting resonantly with the 
cavity field but now we consider the cavity to be lossy, decaying at a rate k. We assume that a detector D is placed outside the 
cavity in a way that it may monitor the cavity decay. 

The description of damping in quantum optical systems is usually described using master equations and its solution gives 
the time evolution of the system when the decay is not observed. However, the time evolution under continuous observation 
of photon counts may be adequately described by a pure state that evolves according to a non-Hermitian Hamiltonian. This 




B. Phonon statistics and continuous observation 
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approach is known as quantum jumps or quantum trajectories 12511 formalism. The idea of continuous observation of decaying 
channels in systems consisting of atoms or ions and cavities has proved itself to be useful to perform legitimate information 
processing tasks as teleportation l23l and maximally entangled state generation IztIE^I or quantum gates l29ll . for instance. We 
saw that the time evolution of the system under Hamiltonian Q, and the realization of a measurement on the electronic state 
may be used to generate the states (|6j. Now, instead of measuring the atomic state, we will show that a measurement of the 
photon outside the cavity collapses the state of the system in a state that keeps much of the characteristics of the state (|6}, namely 
oscillatory behavior of the distribution Pm as well as its narrowing and broadening (l^. For the sake of simplicity, we assume 
that the detector D is perfect. Otherwise we would just have to account for a finite probability that the detector fails in detecting 
an event of leaking of a photon, what would lead us to a description in terms of density matrices rather than state vectors. The 
time evolution of the system conditioned to a no photon decay is given by 



where 



Hcff = —in. h ng 



dt 



??2(l + 2ata) 
2 



(8) 



{a^b^ +a+b). 



(9) 



It is worth noticing that once the Hamiltonian (|9} is not Hermitian, the norm of \ip{t)) is not constant in time. So, it must be 
normalized in order to allow one to correctly evaluate any property of the system. It is clear that if the initial state of the system 
is the same as before, namely, \ip{0)) ~ \a)v\0)c\e), the solution of equation (|8} may be written as 



1^(0) = ^a„i{t)\m,0,e) + b,nit)\m, l,g) 



(10) 



Substituting (I10> and (|9} into (|8} one obtains two coupled differential equations that may be easily solved and the result is given 
by 



amir) = c™(0)e-r^/'^ (^C(t) 
6™(r) = -4ic„(0)e-r-/4 



r 



- 16A- 



2 

LD 



Sir) 



Sir) 



where c,„(0) are the coefficients of the expansion of the initial coherent state in the Fock basis, T = n/g, r ^ gt, and 



C(r) = cosh(./r2-16Ai^r/4) 



Sir) = sinh(./r2-16Ai^r/4) 



(11) 

(12) 
(13) 



Now, we suppose that one photon is detected outside the cavity. This event would correspond to the destruction of one photon 
leading the system to state b\tpir)). Again, we remember that since the time evolution is not unitary the state must be normalized 
after this jump. In our case, the resulting state would be \tpir))d = |$(r))t, |0)c|e), i.e., a disentangled state having a normalized 
motional part given by 



Mr)). = 



bmir) 



m=f) 



|to)„. 



(14) 



Before investigating the statistical properties of that state, it is important to calculate the probability for a photon emission 
because it is related to our probability of success in generating |<I>(t)) u. The probability that at least one jump occurs between 
the initial instant and the subsequent instant r is given by P(t) = 1 — (?/'(t) |?/)(t)), where \4'iT)) is the state in equation (llOt 
with the coefficients M \\ . In figure|3]we have a plot showing the behavior of P(t) using parameters that are close to the ones in 
a current experimental situation, i.e. the weak coupling regime. 

Let us now start the analysis of the statistical properties of the vibrational state |$(r))i,. The ion started in a coherent state 
which has a Poissonian distribution. We can describe its narrowing (or widening) via the normalized variance (also know 
as the Fano factor) defined as = (m^/m) — ifi where m and m? are the first and second moments of the distribution 
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FIG. 3: Probability of detection of one photon outside the cavity. The system parameters are F = 1, -q — 0.05, and a — 2. This probability 
tends to one for higher values of r. 
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FIG. 4: Time evolution of the normalized variance. The system parameters are T — 1, i] = 0.05, and a = 2. 

P,n = |(m|$(r)±)„p, respectively. Values of cr < 1 indicate sub-Poissonian, cr > 1 super-Poissonian, and (7=1 Poissonian 
statistics. The time evolution of (j{t) is shown in figure|4] The original Poissonian distribution naturally evolves to either sub 
or super-Poissonian values. These changes in the width of the distribution could be observed even in a bad cavity that has a 
decay rate k comparable to the coupling constant g, as we can see in figure|3 We would like also to show that strong signatures 
of nonclassical behavior, such as the oscillations in the phonon distribution P,„ at times when the statistics is Poissonian, still 
persist in the weak coupling regime. This may be seen in figure |5] where we show the distribution at a time t = 3.29 and with 
77 = 0.05. Based on those considerations we conclude that general properties of coherent states superpositions, which arise in 
the lossless case, would still persist in our more realistic setup. This means that our proposal could be useful for the experimental 
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FIG. 5: Phonon distribution at r = 3.29 with system parameters F = 1, 77 = 0.05, and a — 2. 
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investigation of certain nonclassical features. 



IV. CONCLUSIONS 

We have investigated several aspects of the dynamics of a trapped ion inside a cavity. Firstly we have considered a situation 
in which the unitary time evolution leads to a global entangled state involving superposed motional coherent states, Fock photon 
states and the two internal electronic states. After the measurement of the internal state of the ion in a specific interaction time, 
the generation of quantum superposition of coherent states of motion of the ion is accomplished. Two different states may be 
generated (either or |$-)) depending on the result of the measurement of the internal ionic state. The main requirement 
for such generation is the strong coupling regime where the system may perform Rabi oscillations in the lifetime of the cavity 
photon. In the second part of our paper we consider the influence of cavity decay in the ionic dynamics. In fact that represents 
itself a generation method, since a nonclassical state results from the dissipative evolution even with a photon decay rate of 
the same order as the ion-cavity coupling (weak coupling regime). The cavity is continuously monitored by a detector, what 
causes the state of the system to be pure at any time. The measurement of the internal electronic state in the former suggestion 
is replaced now by the counting of a photon leaking out of the cavity. This collapses the entangled global state of the system 
onto a product state. Even though the cavity is not ideal, the ionic motional state still retains (after the photon decay) important 
nonclassical features that characterize quantum superposition of coherent states, such as, for instance, changes in the variance of 
the phonon distribuion (sub or super-Poissonian statistics) as well as its oscillatory behavior. 
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